Draft version November 6, 2012 

Preprint typeset using I^T^X style cmulatcapj v. 5/2/11 



A SEARCH FOR TRANSITS OF GJ 581e AND CHARACTERIZATION OF THE HOST STAR VARIABILITY 

USING MOST* SPACE TELESCOPE PHOTOMETRY 

Diana Dragomir 1 , Jaymie M. Matthews 1 , Rainer Kuschnig 2 , Jason F. Rowe 3 , Brett J. Gladman 1 , David B. 
Guenther 4 , Anthony F. J. Moffat 5 , Slavek M. Rucinski 6 , Dimitar Sasselov", Werner W. Weiss 2 

Draft version November 6, 2012 

ABSTRACT 

The GJ 581 system has been amply studied since its discovery in 2005: the number of known planets 
in the system has increased and their orbital parameters are among the most precisely determined for 
radial velocity detected exoplanets. We have acquired MOST space-based photometry during 2007 
and 2009, with the aims of measuring the stellar variability and searching for transits of GJ 581e, 
respectively. We quantify our sensitivity to shallow transit signals using Monte Carlo simulations, 
and perform a transit search within the 3<r transit windows corresponding to both the circular and 
Keplerian orbit ephemerides. Our analysis rules out transits for a planet with an orbital period of 3.15 
days (GJ 581 e) having a radius larger than 1.62 (or a density lower than 2.39 g cm -3 for an orbital 
inclination of 90°) to 2a confidence. Thus, if the planet transits, we can exclude hydrogen, helium and 
water theoretical model compositions. The MOST photometry also allows us to rule out transits of GJ 
581b within the Keplerian orbit-derived transit window for impact parameter values smaller than ^0.4 
and confirm previous results which exclude transits for this planet within the circular orbit-derived 
transit window, for all plausible interior compositions. We find that the stellar brightness of GJ 581 
is stable to within 1%, a characteristic which is favourable to the development of life in the habitable 
zone of the system. In the 2009 photometry, we detect a stellar signal with a period of 5.586 ± 0.051 
days, which is close to the orbital period of GJ 581b (P =5.37 days). However, further monitoring of 
the system is necessary to verify the nature of this variation. 

Subject headings: planetary systems - techniques: photometric - stars: individual (GJ 581) 



1. INTRODUCTION 

Radial velocity (RV) and transit searches are discov- 
ering progressi vely less massive an d small er exoplanets. 
The HARPS ( IMayor et aT1[2003| , Keck ( IButler et a~ 
1996 1 several other RV surveys have brought forth dozens 
of planets wit h minimum masses < 10 Mg, while the Ke- 
pler mission ( |Borucki et al.|p010[ ) has so far yielded a 
wealth of super-Earth candidates with radii <4 R®. Fur- 
thermore, Kepler has already detected three planetary 
candidates with radii smaller than that of the Earth in 



a single system ( Muirhead et al 



2012). 



There are currently a handful of super-Earths for which 
both the masses and radii have been measured. By deter- 
mining the density of these planets, of which none exist 
in the Solar system, it becomes possible to constrain their 
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composition and eventually their formation. Thanks to 
their transiting nature it also becomes possible to per- 
form spectroscopic studies of these planets. Such stud- 
ies, rendered possible by the large transit depth which is 
induced by the small size of the host star relative to that 
of the tra nsiting planet, now exist for the super-Earth 
GJ 1214b flBean et al.|2010||Bean et al.|2011[|Croll et al 
Bertaet al.||20l2| ). " 

58 1 , a small M2 . 5 dwarf (i?*=0.299 R^;|von Braun 



2011 



G. 



et al 



et al 



( 2011 )), is known to host up to 6 planets 
20Q'l[l dry e~ 



al.|2010p . Ho 

t and g nas b< 



Bbnfils 

Mayor et al.||2009||Vogt et 



wever, we note that the existence of planets 
gTias been disputed and arguably ruled out by sev- 



eral independent a nalyses ( Gregory[2011 Forveille et al. 
201 lj Tuomi||2011| ). This system has been observed and 



studied at length. The published HARPS radia l velocity 



measureme nts span approximately seven years ( Forveille 
et al. 2011 ), and the HIRES observations cover more than 



a decade. The orbital param eters of four planets i n the 



system are well constrained (Forveille et al. 2011 ). At 
least three of those planets (c, d and e) have minimum 
masses in the super-Earth regime. The habitability of 
c and d has been assessed by several groups, and while 
the former orbits too close to the star to sustain liquid 



water on its surface (von Bloh et al. 2007 Selsis et al 



20071, GJ 581 d was found to reside in the habitable zone 



of the system ( von Bloh et al.||2007[ |Wor dsworth et al. 
20TTl|barnes et al.|2009[|von Paris et al.|2011[|Kalteneg 



ge r et al.||201ip . 

1'he GJ 581 system was observed by th e MOST (Mi 



crovariability and Oscillatio ns of STars; | Walker et al.| 
2003 Matthews et al. 2004) space telescope in order to 
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TABLE 1 

Published O rbital Parameters for GJ 581e 

FROM [FpRVEILLE ET AL~] ( |20 11 [ ) 



TABLE 2 

Published O rbital Parameters for GJ 581b 

FROM [FPRVEILLE ET Al!] ( |201 1 [ ) 



Parameter 


Circular model 


Keplerian model 


Parameter 


Circular model 


Keplerian model 


P (days) 


3.14941 ± 0.00022 


3.14945 ±0.00017 


P (days) 


5.36864 ± 0.00009 


5.36865 ± 0.00009 


To a (HJD) 


2454749.026 ± 0.056 


2454750.31 ± 0.13 


T a (HJD) 


2454751.536 ±0.012 


2454753.95 ± 0.39 


e 




0.32 ±0.09 


e 




0.031 ± 0.014 


w (deg) 




236 ± 17 


u (deg) 




251 ± 26 


K (ms" 1 ) 


1.754 ±0.180 


1.96 ±0.20 


K (ms" 1 ) 


12.72 ±0.18 


12.65 ±0.18 


M p sin i (M®) 


1.84 


1.95 


M p sini (M®) 


15.96 


15.86 


a (AU) 


0.028 


0.028 


a (AU) 


0.041 


0.041 



a Predicted mid-transit time. 

determine its variability and search for transits of GJ 
581e. These observations were acquired as part of a pro- 
gram aiming to search for transits of known RV super- 
Earth candidates using MOST. The telescope is opti- 
mized for bright star (V < 10) photometry, a characteris- 
tic it shares with radial velocity surveys. Even though GJ 
581 is slightly fainter (V = 10.6) than the usual MOST 
targets, the opportunity to obtain continuous coverage of 
more than one orbital cycle of the planet (thus eliminat- 
ing the chance that a transit might be missed) motivated 
the observations. 

GJ 581e, the innermost planet in the system, orbits its 
host star with a period of 3.15 days and has a a priori 
transit probability of 5%. This value corresponds to the 
geometric probability that the planet crosses the disk of 
the star and in the case of a circular orbit it is derived 
from the planet's semi-major axis and the stellar radius. 
GJ 581e, whose published properties can be found in 
Table 1, is the planet with the smallest radial velocity- 
derived minimum mass to date. If the planet is found 
to transit, its size can be determined and spectroscopic 
studies of its atmosphere become possible. 

The duration and nearly contiguous nature of the 
MOST observations presented in this paper also invite 
a search for transits of GJ 581b (P = 5.37 days), re- 
gardless of the outcome of our GJ 581e transit search. 
Indeed, even if GJ 581e does not transit, GJ 581b may 
still transit if their orbits are not co-planar. Before the 
existen ce of the other three planet s in the system was 

(2006) had ruled out tran- 



known, Lopez-Morales et al. 

sits of GJ 581b using an orbital solution based on the 
assumption of a circular orbit. However, when the three 
subsequently discovered planets are taken into account 
and their orbital eccentricities are allowed to vary dur- 
ing the fitting process, the resulting solution shifts the 
predicted mid-transit t ime of GJ 581b outside the ob- 
servations acquire d by Lopez-Morales et al. (20061, as 
was determined by Forveille et al. (2011) 



The published 
o rbital parameter s of G J 581b are found in Table 2. 



Forveille et al. 



(20111 performed a transit search for 
GJ 581e using photometry from the 2.5 m Isaac Newton 
Telescope that was obtained during one transit window 
between -1 and 2.3 a from the predicted mid-transit time 
(a total of 6.5 hours) for a circular orbit. They have found 
no evidence of a transit. In this paper, we present two 
sets of nearly contiguous MOST space-based photometry 
of GJ 581, obtained in 2007 and 2009 (section 2). The 
2009 data set covers almost four full orbital cycles of GJ 
581 e (section 3.1) and thus allows us to perform a more 
complete transit search using the ephemeris for the cir- 



a Predicted mid-transit time. 

cular orbital solution. In addition, we are able to search 
for transits during the predicted transit window^] corre- 
sponding to the eccentric orbital model. We quantify our 
sensitivity to shallow transits via Monte Carlo tests (sec- 
tion 3.2). In sections 3.3 and 3.4 we place limits on the 
size and mass of GJ 581e and discuss the implications. 
We use the MOST photometry to search for transits of 
GJ 581b within its circular and Keplerian orbit-predicted 
transit windows in section 4. In section 5, we briefly ad- 
dress the cases of GJ 581c and GJ 581d. Finally, we 
combine this data set with a longer light curve obtained 
in 2007 to carry out a stellar variability analysis for GJ 
581 (section 6). We conclude in section 7. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained photometry of Gliese 581 from the MOST 
satellite during 2007 April 26 to May 24, and during 2009 
May 8-20. MOST is a micro-satellite which carries a 15 
cm optical telescope feeding a CCD photometer. It is 
in a Sun-synchronous polar orbit with a period of 101.4 
minutes, which allows it to monitor stars in a Continu- 
ous Viewing Zone (CVZ) without interruption for up to 8 
weeks. The CVZ covers a declination range of -18° < 5 < 
+36°. Stars brighter than V ~ 6 are observed using a 
Fabry microlens to project onto the CCD an image of the 
telescope pupil illuminated by the target. Fainter stars 
are observed in Direct Imaging mode, in which the defo- 
cused images of the stars are projected onto the CCD. 

In 2007, MOST observed GJ 581 for 50 days (Figure 
1), but alternated between this star and other targets 
during this period. Two stretches of time during which 
the sampling cadence was significantly higher than for 
most of the run can be seen in Figure 1. The intention 
of these two more concentrated "windows" of coverage 
was to search for transits of the super-Earth GJ 581c. 
Unfortunately, due to a mistranslation of the radial ve- 
locity ephemeris, these high-cadence sections of the 2007 
light curve do not overlap with predicted transit windows 
for this planet, nor with any predicted mid-transit times 
for the other known planets in the system. For this rea- 
son, and because the sampling outside of these windows 
is not sufficiently contiguous, the 2007 photometry was 
not effective for a transit search. Instead, it was com- 
bined with the 2009 data to assess the variability of the 
host star (see section 6). In 2009, GJ 581 was observed 
contiguously during approximately half of every MOST 

9 The transit window is the time span during which a transit is 
predicted to occur, calculated from the uncertainties on the orbital 
period and those on the predicted mid-transit time. 
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Fig. 1. — Reduced MOST photometry acquired during 2007 April 25 to June 14. The tire-track pattern of the data is due to the fact 
that MOST alternated between GJ 581 and other stars during these observations. Sec section 2 of the text for details. 
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Fig. 2. — Top: Reduced MOST photometry acquired during 2009 May 8-20. Bottom: Same, but with the modulation visible in the top 
panel removed (see section 2 for details). The tire-track pattern of the data is due to the fact that MOST alternated between GJ 581 and 
another star during these observations. See section 2 of the text for details. 



orbit (54 out of every 101.4 minutes), for nearly 12 days 
(Figure 2). The alternating target was monitored during 
the remaining 47.4 minutes. Both the 2007 and 2009 ob- 
servations were acquired in Direct Imaging mode, with 
individual exposure times of 3.0 s. Sets of 8 individual 
frames were stacked on board the satellite before being 



downloaded to Earth, leading to a total integration time 
of 24 s per data point. There is also a modest readout 
overhead of 3.2 seconds, resulting in a sampling rate of 
27.2 s (outside of the interruptions for the alternating 
target). 

The raw data were reduced using aperture photometry. 
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The re duction pipeline (described in detail in Rowe et al 
20081 corrects for cosmic ray hits and stray light from 



scattered Earthshine, which varies with the period of the 
satellite. The 1-day period variation (due to the Sun- 
synchronous orbit of MOST) with an amplitude of 1.3 
mmag was also filtered from the data. 

In addition, for the 2009 data set three short sections 
of the light curve (at around 3418.2, 3420.5 and 3424.6 
days in Figure 2) affected by high levels of stray light 
were conservatively excised. After these steps, a modu- 
lation remains in the light curve (Figure 2, top panel). A 
possible origin of this modulation is discussed in section 
6. This trend was modelled with a cubic spline inter- 
polated from the data binned every 500 minutes, and 
removed from the light curve before the transit search 
and transit-injection tests. The bin size was chosen to 
be large enough not to interfere with any transits of GJ 
581 e (predicted edge-on transit duration = 1.2 hours) 
and small enough to remove the slowly varying trend. 
The final 2007 and 2009 time series contain 7667 and 
17312 points, respectively. 

3. GJ581e 

3.1. Transit Ephemeris and Predicted Characteristics 



Our approa ch is inspired by the m ethod s used in|Croll 



et al. 



We use the ephemerides published in |Forveille et al. 
(2011) to search for transits of GJ 581e. The authors 



derived best-fitting orbital models for planets b, c, d and 
e in two ways: by assuming circular orbits (where the 
eccentricity is fixed to 0.0 for each planet) and by fitting 
fully Keplerian orbits (where the eccentricity is allowed 
to float). While the values of the orbital period for GJ 
581e in the two cases are nearly identical, an eccentricity 
of 0.32 ± 0.09 is found for th e Keplerian fit (see Table 
1 for the |Forveille et al. |2011| orbital parameters of GJ 
581e). In addition, Mp sin i is slightly larger for the Ke- 
plerian orbit fit (1.95 M s ) than for the circular orbit 
model (1.84 M s ), though no uncertainties are reported 
for these values. The predicted mid-transit times (TO) 
and associated uncertainties also differ significantly be- 
tween the two cases, leading to different phase locations 
and sizes for the transit windows. The coverage of the 
MOST photometry is adequate for monitoring the tran- 
sit windows corresponding to both solutions. 

For the circular orbit transit search, the photometry 
were folded on TO =2454969.489. The la uncertainty 
on this value is 0.071 days, so the size of the 3<r tran- 
sit window is 11.48 hours (including a predicted tran- 
sit duration of 1.2 hours). In the Keplerian orbit case, 
TO =2454970.77 was used to fold the time series. The as- 
sociated lcr uncertainty is 0.14 days, leading to a 3<r tran- 
sit window of 22.03 hours (including a predicted transit 
duration of 1.6 hours). 

We estimate radii for GJ 581e for each of three theoret- 
ical compositions: hydrogen, water ic e and silicate. We 
use the models of Seager et al. ( 2007 1 and the planetary 



mass obtained by assuming an impact parameter (b) of 
0.0 (i = 90°). The derived radii are 3.70, 1.68 and 1.25 
J?^, for a homogeneous hydrogen, water ice and silicate 
planet, respectively. 

3.2. Monte Carlo simulations 

Throughout this section we only use the 2009 data set 
because the sampling cadence of the 2007 photometry is 
too low for a transit search of a super-Earth planet. 



(|2007a|), iCroll et all (p007bl and IBailard et al 



120101. However, our analysis differs in that we search 



for transits of a RV-detected planet with a well-known 
period rather than additional planets in known transiting 
systems. Indee d, the orbital period o f GJ 581 e is known 
to within 19 s ( Forveille et al.|[20lT |. As a consequence, 
for our Monte Carlo simulations and transit search we 
vary the planetary radius (R p ), orbital inclination (?) and 
orbital phase at which the transit may occur (0 < </> < 1), 
but not the orbital period. Carrying out simulations for 
values of i other than 90° allows us to quantify the effect 
of the sampling cadence of the MOST observations on 
our sensitivity to transits of various durations. 

To quantify our detection limits for the Gliese 581 pho- 
tometry, we injected simulated limb-darke ned transits 
(using the models of |Mandel fe Agol|[2002| in the 2009 
light curve and attempted to recover them by fitting a 
box-shaped transit model and minimizing \ 2 ■ For the 
simulated transits, we interpolated limb darkening coef- 
ficients (ci =1.600, c 2 =-1.491, c 3 =1.209, c 4 =-0.398) 
from tables generated for the MOST bandpass (A. Prsa, 
pri vate communication) and based on the Kurucz mod- 
els ( |Castelli fc Kurucz|2004| ), assuming a stellar effective 
temperature ('Je//) of 3500 K, surface gravity (log g) of 
5.0 and metallicity ([Fe/H]) of -0.15. The actual values 
for GJ 581 are T eff = 3498 ± 56 K, log q = 4 .96 ± 0.08 
and [Fe/H] = -0.135 ( |von Braun et al.l|2011[ ). We used 
a grid evenly spaced in radius [Rp-i n j) and orbital in- 
clination (ii n j). 



The grid consists of 30 i?p_i„, values 



ranging from 1 to 4 i?© and 4 ii n j values ranging from 
87.5 to 90° (0.0 < b < 0.9). For each combination of ra- 
dius and inclination, we insert transits at 100 randomly 
distributed phases {4>inj) spanning the 3<r transit win- 
dow for the circular ephemeris. Thus we search a total 
of 12000 light curves. 

We sample the x 2 space at 90 values of the planetary 
radius (equally spaced from 0.85 to 4.5 i?© ), 6 values 
of the orbital inclination (equally spaced from 87.2 to 
90°) and between 130 and 2000 values of the orbital 
phase (equally spaced within the 3er transit window). For 
each model, we compute the change in \ 2 relative to the 
constant flux hypothesis. The model with the largest 
value of A X 2 % = 100%( X 2 - xD/xl (where X t and X l 
correspond to the transit and constant flux model, re- 
spectively) becomes the best-fit solution for a particular 
light curve. We call the values of R p , i and (j> associated 
with the best-fit solution R rec , i rec and </> rec , respectively. 
Since the goal is to identify a statistically significant tran- 
sit rather than to precisely determine its parameters, we 
consider a transit as successfully detected if the following 
two conditions are satisfied: | (R r ec/Rinj)\ — 1 < 20% and 
< 0.01 (corresponding to ~ 45 min). Rj,. 



are the planetary radius and orbital phase of the 



|0rec 

and i 

injected transit signal. The results of our Monte Carlo 
simulation indicate that we do not always accurately re- 
cover the impact parameter of the injected signal. This 
is not surprising for such shallow and short (1.2 hours 
for 6 = 0) transits. As a consequence, we do not use an 
orbital inclination condition as a criterion for detection. 

For b = (i = 90°), we recover the transits correspond- 
ing to a planetary radius >1.42 R® at least 68% of the 
time. We recover injected transits for R p > \.Q2R^ and 



5 



for R p > 1.93i?® at least 95% and 99% of the time, re- 
spectively. For larger values of b these confidence limits 
tend to shift to larger radii, though only mildly so for b < 
0.6. This is because the transit duration is shorter and 
there are less in-transit data points, making the signal 
more difficult to detect. Further, for larger values of b 
the planet would cross nearer the limb of the star, thus 
blocking out less light and leading to a shallower transit 
than in the case of b = 0. In light of these two argu- 
ments, we note that while transits with Rp < 4i?^ and 
b = 0.9 do not meet our detection criteria with sufficient 
statistical confidence, deeper transits do. Our detection 
limits as a function of b and R p are shown in Figure 3, in 
the form of 68%, 95% and 99% confidence contours. A 
value of b = corresponds to an edge-on transit, while 
6 = 1 represents a grazing transit. 
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Fig. 3. — Detection limits for injected transits as a function of 
planetary radius (Rp—i n j) and impact parameter for GJ 

581e. The solid points correspond to the values of Rp_ in j and 
bi n j at which we tested our sensitivity. The solid, dashed and 
dotted lines indicate the 68%, 95% and 99% confidence contours, 
respectively. 



3.3. Transit Search 

We perform the transit search described above on the 
2009 light curve (within the 3<r transit window corre- 
sponding to the circular orbit solution). We do not find 
any transit candidates for a planet with R p >1.42 
which is our la detection limit (for b — 0). 

Our best transit candidate (for which A\ 2 % = 0.54% 
was largest) has the following properties: R p = O.912i?0, 
b = 0.40 and cj) = 0.058, and therefore it occurs near the 
positive edge of the 3c window (the predicted mid-transit 
time corresponds to <j) = 0.0). This signal is significantly 



below our detection limit and over five times smaller than 
the rms of the phased data within the 3cr transit window. 
Such a planet would have a density (p) of 13.4 g cm -3 , an 
unrealistic value as it is greater than the maximum iron 
fraction limit for a rocky planet ( Marcus et al.|20i0 ). For 
these reasons, we consider this candidate very unlikely 
to be a transit signal. The phased observations spanning 
the 3er transit window and the best-fit model for this 
candidate are shown in Figure 4 



Using the HIRES-based results from |Vogt et al | ( |2010| ), 
Andrae et al. ( |2010 ) find a possible discrepancy between 
the residuals of a circular orbit solution for the four plan- 
ets, and the expected Gaussian distributed residuals if 
the best-fit model was the true model. They suggest 
this discrepancy may indicate that the planetary orbits 
are eccentric. Further, when fitting a Kepl erian model 
to the HARPS RV data (see section 3.1), Forveille et 



al. (2011) obtain significant values for the eccentricity 
01 planets d and e, and different mid-transit times for 
all planets than in the circular model case. Motivated 
by these result, we also carry out a transit search for 
GJ 581e within the 2>a transit window computed from 
the orbital parameters associated with the Keplerian 
model. The rms=0. 00425 (in units of relative flux) of 
the phased photometry within this window is marginally 
smaller than that of the circular solution transit window 
(rms=0. 00431). Thus our sensitivity to transits should 
be very similar or marginally better. We find no transit 
candidate near or above the detection limits reported in 
Figure 3. 

The best candidate (Ax 2 % = 0.92%) in this case 
represents a planet with R p = 1.19i?®, b = 0.79 and 
<f> = —0.006, very close t o the predicted mid-t ransit time. 
With a mass of 1.95 M e ( TForveilfe et al.|2011[ ), this planet 
would have p =6.41 g cm -3 which would make it slightly 
denser than 55 Cnc e and suggest a silicate composition 
model (Seager et al. 20071. Nevertheless, the statistical 
significance of this signal is too low relative to our de- 
tection limits to make it a likely candidate. The phased 
observations showing this candidate can be seen in Fig- 
ure 5. 

3.4. Constraints and Implications 

Even when a super-Earth transits and its radius can 
be measured, it is still difficult to gain insight into its 
detailed composition. Nonetheless, useful limits can be 
placed on its interior structure when both it s mass and 

Unfortu- 



2010) 



ransit signal for 



radius are available (Rogers & Seager 
nately, we do not detect a significant 
GJ 581e. If the planet transits, then we can use the 
limits determined in section 3.2 to constrain its size as 
follows. 

To la confidence, we exclude pure hydrogen, hydro- 
gen/helium, pure water ice and 50% water ice composi- 
tions for detectable (b < 0.6) transiting orbital configu- 
rations. To 2a confidence, we rule out compositions less 
dense than 75% water for b < 0.6. We exclude hydrogen 
and hydrogen/helium compositions with 3er confidence. 
Quantitatively, for Mp — 1.84M® (based on the circular 
orbit solution), we rule out densities lower than 3.54 g 
cm~ 3 , 2.39 g cm -3 and 1.41 g cm~ 3 with la, 2a and 3a 
confidence, respectively. 

There are so far no transiting planets with a measured 
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Fig. 4. — The best candidate transit signal for the circular model 
ephemeris of GJ 581e. Phase 0.0 corresponds to the predicted 
mid-transit time. Top: Unbinned photometry folded at the orbital 
period of 3.15 days. The orbital phase range of the plot is equiv- 
alent to the 3<r transit window. Bottom: Zoom-in of the folded 
photometry, binned every 0.001 orbital phase. In both panels, the 
best-fit box-shaped transit model is shown in red. 
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Fig. 5. — The best candidate transit signal for the Keplerian 
model ephemeris of GJ 581e. Phase 0.0 corresponds to the pre- 
dicted mid-transit time. Top: Unbinned photometry folded at the 
orbital period of 3.15 days. The orbital phase range of the plot 
is equivalent to the 3<r transit window. Bottom: Zoom-in of the 
folded photometry, binned every 0.001 orbital phase. In both pan- 
els, the best-fit box-shaped transit model is shown in red. 



mass smaller than 2 M®. Kepler-llf comes close, at 2.3 
M ffi (determined from dynami cal fits of the Kepl er-11 
system) and a radius of 2.6 i?® (Lissauer et al. 2011 ). GJ 
1214b orbits a star much closer in spectral type (M4.5) 
to GJ 581, but has a mass of 6.4 M m and a radius of 



2.7 i?® (Charbonneau et al. 20091. These two super 



Earths a re relatively low-d ensity planets, wit h envelopes 



of water ( Berta et al.|2012 ) and/or hydrogen ( [Lissauer et 
al. 2011). If GJ 581e transits and has a similar density, 
then we would have detected the transit in the MOST 
photometry. However, if a transiting GJ 581e is more 



=5.9 g cm 

-3 



-3. 



Winn 



similar to t he denser 55 Cnc e (p 
et al.| |201l|), Kepler-lOb (p =8.8 g cm 
2011|), or Corot-7b (p =10. 4 g cm" 3 , assuming a mass of 



Batalha et al. 



Hatzes et aL|2011[ ), then we would not have de- 



7^M 

tected its transit using our data. Nevertheless, 55 Cnc e, 
Kepler-lOb and Corot-7b receive significantly more ra- 
diation due to both their proximity to their host stars 
and the higher temperatures of those stars, compared to 
Kepler-llf, GJ 1214b and GJ 581e. Therefore, we expect 
any atmospheres or water envelopes these highly irradi- 



ated planets may have acquired to be mostly or entirely 
evaporated, thus possibly explaining their higher densi- 
ties. 



Mayor et al. ( 2009 ) have performed a stability analysis 
of the <JJ 581 system in order to determine constraints 
on its orbital inclination. Assuming co-planar orbits and 
allowing for non-zero eccentricity, they found that GJ 
581e was ejected for inclinations less than i — 40°. This 
lower bound leads to a maximum mass for GJ 58 le of 
2.86 M® based on the Mp sin i value from the circular 
solution, and 3.03 -M®, based on the Keplerian solution. 
A co-planar orbital configuration with i > 40° , even if it 
is not edge-on, should still give rise to observable thermal 
phase variations of the order of at least 10 ppm in the 
infrared, dominated by the thermal radiation of planets 
b (see section 4) and e. A signal of this am plitude is 
lik ely wi thin reach of NIRCam on the JWST ( Seager et 
aLl [2009| ). 



4. GJ 581b 

4.1. Transit Ephemeris and Predicted Characteristics 
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We use the Forveille et al. (2011) ephemerides for GJ 
581b to determine its 3a transit window for both the 
circular and Keplerian orbit cases, in the same man- 
ner we did for GJ 581e. Although the authors find a 
nearly circular orbit (e = 0.031 ± 0.014) for GJ 581b it- 
self in the Keplerian orbit case, the planet's predicted 
mid-transit time is nevertheless affected (and shifted rel- 
ative to the circular case) due to the overall best-fitting 
solution which does assign significantly greater than 
eccentricities to planets e and d. The minimum mass 
of GJ 581b is not very different in the circular case 
(Mp sin i — 15.96M®) compared to the Keplerian case 
(Mp sin i = 15.86M®), and is very close to that of Nep- 
tune (M P = 17. 15Mffi) 



We find that Lopez- Morales et al. (2006) have con 



clusively ruled out transits of GJ 581b within its ~ 2a 
transit window for the circular orbit ephemeris. They ex- 
clude planetary radii greater than 1.5R^, a value which, 
given the minimum mass of the planet, would suggest 
a highly implausible object composed entirely of iron. 
However, to our knowledge the transit window arising 
from the Keplerian fit has not yet been searched. 

For completeness (the MOST observations allow the 
search to be carried ou t beyond the 2a window) and 
to confirm the results of Lopez-Morales et al. (2006), we 
perform this analysis for the circular orbit window as 
well, using TO = 2454966.303 to fold the photometry. 
The la uncertainty on this value is 0.016 days. The 
3a transit window (including a predicted edge-on transit 
duration of 1.4 hours) is thus 3.7 hours long. For the 
Keplerian orbit window we use TO = 2454963.34 to fold 
the data, and the associated la uncertainty of 0.39 days 
leads to a transit window length of 56.2 hours (or 2.34 
days). 

4.2. Transit Search and Constraints 

The folded photometry within the 3a transit window 
for the circular and Keplerian orbit cases are shown 
in Figures 6 and 7, respectively. We considered and 
overplotted transit models corresponding to three pos- 
sible planetary radius values. The dashed line repre- 
sents the minimum radius for a planet of mass equal 
to the minimum mass of GJ 581b with the maximum 
iron core mass fraction attainable from collision-induced 
mantle strip ping for a planet (Rp = 1.9R, 



Marcus 



et al. 2010). A more realistic transit depth is indi 
cated by the dotted line, which corresponds to a silicate 
planet (Rp — 2.3i?®). Finally, the solid line represents 
a planet with a density equivalent to that of Neptune 
(Rp = 3.8i?e). 

The phase coverage for GJ 581b is incomplete even 
when the ^2 orbital cycles spanned by the observations 
are phase-folded because of the sampling cadence of the 
MOST photometry of this system (see section 2). This 
sets limits on the duration of a transit detectable with the 
MOST data, depending on where it would occur within 
the window. 

Thus, by visual inspection of the phased light curve 
within the predicted transit window, we can rule out 
transits corresponding to radii associated with any plau- 
sible composition (Rp > 2.3i?®) for GJ 581b, but only 
for impact parameter b < 0.4 (corresponding to a tran- 
sit duration greater than ~ 90% of that of an edge-on 
transit). These limits apply to both the circular and Ke- 
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Fig. 6.— The 2009 photometry, folded at the period of GJ 581b 
(5.37 days). The orbital phase range of both top and bottom plots 
is equivalent to the 3<r transit window based on the circular or- 
bit ephemeris. Phase 0.0 corresponds to the predicted mid-transit 
time. Box-shaped predicted transit signatures for three planetary 
radius values are overplotted in blue. The dashed, dotted and 
solid lines represent Rp = 1.9,2.3 and 3.8R^, respectively. Top: 
Unbinncd photometry. Bottom: Vertical zoom-in of the folded 
photometry, binned every 0.0006 orbital phase. 



plerian orbit cases. It is difficult to set more precise (and 
applicable to the entire transit window) limits on the 
duration of a detectable transit because the sampling ca- 
dence and phase coverage of the phased photometry are 
uneven throughout the transit window. 
For the circular orbit window - based on the results of 

we conclude 



Lopez-Morales et al. 
that transits are un 



(2006|) and our own 
«5ly 



o occur. For the Keplerian 
case, given our observations, transits may still occur for 
0.4<6<1. 

The lower limits on the mass of GJ 581b thus do not 
chang e . The dynamical simulations performed by |Mayor] 
(2009) (under the assumption of stable and co- 



et al. 



planar orbits for planets a to d) lead to an upper limit 
of 24.83 M® for the circular orbit case and 24.67 M® for 
the Keplerian orbit case. 

5. GJ 581c AND GJ 581d 

The third and fourth known planets in the GJ 581 
system have orbital periods of 12.92 and 66.6 days, re- 
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Fig. 7. — The 2009 photometry, folded at the period of GJ 581b (5.37 days). The orbital phase range of both top and bottom plots 
is equivalent to the 3cr transit window based on the Keplerian orbit ephemeris. Phase 0.0 corresponds to the predicted mid-transit time. 
Box-shaped predicted transit signatures for three planetary radius values are overplotted in blue. The dashed, dotted and solid lines 
represent Rp = 1.9,2.3 and 3.8R&, respectively. Top: Unbinned photometry. Bottom: Vertical zoom-in of the folded photometry, binned 
every 0.0006 orbital phase. 



spectively. They are both super-Earth candidates, with a 
minimum mass of 5.4 M® for GJ 581c and 5.3 Af® for GJ 
581d, in the circular orbit case. Their long orbital periods 
result in low transit probabilities, but having photometry 
spanning several days to several weeks motivated us to 
consider whether the observations may include the pre- 
dicted mid-transit times of these two planets. 

Unfortunately, no GJ 581 d predicted transit times oc- 
cur within the 2007 or 2009 data sets, for either the Kep- 
lerian or circular orbit ephemerides. No GJ 581c transit 
times are predicted within the high-cadence sections of 
the 2007 light curve (visible in Figure 1), but the 2009 
data set does include one value of TO for each of the two 
orbital fits. For the circular case, TO = 2454970.93 with 
a Ict uncertainty of 0.15 days. For the Keplerian case, 
TO = 2454969.7 with a la uncertainty of 1.6 days. 

The main reason we do not carry out a transit search 
beyond visual inspection of the light curve is that the 
MOST photometry only covers one transit window for 
GJ 581c. Most of the possible radius values for this 
planet are too small to produce a transit signal that 
would be believable with a single detection given the 
precision of our photometry. The relatively long transit 
windows, especially for the Keplerian orbit case, would 
further decrease the significance of such a signal. We do 



note that no transit feature with depth larger than the 
rms of the light curve (~ 0.005) is visible in the photom- 
etry. This limit corresponds to a radius of 2.3 i?® in the 
case of an edge-on transit. Given the relatively low mass 
of GJ 581c (5.4 M®) this is not a particularly constrain- 
ing limit on the radius of this planet. If GJ 581ctransits 
the host star, this limit still allows a range of densities 
which include water/helium/hydrogcn planetary compo- 
sitions. 

6. STELLAR VARIABILITY 

The brightness of GJ 581 is stable to less than 1% 
over the four weeks during which MOST observed the 
star in 2007. The 11.5-day data set obtained in 2009 
supports this conclusion. This level of stability agrees 
with the finding that the X-ray brightness of GJ 581 
lies below ROSAT's detection threshold and supports an 
advanced age (> 7 Gyr) for the star, as suggested by 



Selsis et al. ([2007). Thus, it constitutes a factor in favour 
ot tfie potential development of life on any planet residing 
within the habitable zone of the system, and all ows a 



sufficiently long period of time for its occurrence (|Scalo 
eFal][2007| 



We searched the MOST photometry for eviden ce of 
variations indicative of the stellar rotation period. |Vogt| 
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et al. (2010) report a value of 94.2 days, based on 
AFT photometric observations acquired with the 0.36 m 
Tennessee State University telescope and spanning 530 
nights. They find a semi-amplitude of 3 mmag (0.0028 
in units of flux) for this variation. Though it covers only 
31% of the reported stellar rotation period, the MOST 
2007 data is contemporaneous with a subset of the 0.36 
m photometric observations ranging from 2454217.3 to 
2454244.5 (rotational phase 0.84 to 1.1 3 acco rding to 
section 4 and Figure 1 of |Vogt et al] ( |2010| )). The 
MOST 2007 photometry shows a slight long-term cur- 
vature, with a minimum at approximately HJD 2454236. 
Th is cor responds to phase 1.05 in Figure 1 of |Vogt et] 



0.0020 



al. 2010 



In their figure, the minimum of their folded 
photometry can also be observed near the same phase. 
However, our data are insufficient to definitively establish 
whether they support a rotation period of 94.2 days. 

In the 2009 data set we observe a modulation in the 
light curve with a period of about 5 days and semi- 
amplitude of 1.5 mmag (or 0.14% in flux units; see Fig- 
ure 2), which was removed prior to the transit injection 
tests and transit search as described in section 2. The 
amplitude spectrum shows a strong peak at a frequency 
of 0.1790 ± 0.0016 cycles/day (period of 5.586 ± 0.051 
days), as can be seen in Figure 9. This variation is es- 
pecially interesting as it has a period comparable to the 
orbital period of GJ 581b (P = 5.36865 ± 0.00009 days). 
We generate an amplitude spectrum of the 2007 photom- 
etry (see Figure 8) to verify whether the same signal is 
present. A peak of much lower significance is visible at 
a frequency corresponding to a period of 5.365 ± 0.074 
days. In Figure 10, we show the 2007 and 2009 photome- 
try phased to the planet's orbital period, separately and 
together. While the 2007 data is not persuasive, the 2009 
photometry shows a coherent signal. This signal persists 
when the two data sets are combined, and while this is at 
least in part due to the greater number of observations 
in the 2009 photometry which dominate the combined 
data, its persistence is arguably encouraging. 
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Fig. 8. — Amplitude spectrum (in units of relative flux) for the 
2007 MOST photometry. The red line indicates the frequency cor- 
responding to the orbital period of GJ 581b (P = 5.37). 



Regardless of whether this variation persists beyond 
the duration of our 2009 observations, its period is too 




Fig. 9. — Amplitude spectrum (in units of relative flux) for the 
2009 MOST photometry. The red line indicates the frequency cor- 
responding to the orbital period of GJ 581b (P = 5.37). 



short to corre spond to the stella r rotation period, as 



pointed out by Mayor et al 
proximate match between t 



( 2009 ) . Considering the ap- 
le orbital period of GJ 581b 



and the period of the variation, we tested the possibility 
that the variation is due to scattered light from the sur- 
face of the planet. Assuming a conservatively large radius 
of 0.7 Rj U p for GJ 581b, (Rp/a) 2 gives a signal with a 
semi-amplitude of approximately 0.006% or 0.07 mmag. 
This is much smaller than the semi-amplitude of the vari- 
ation we observe. Thus, the variation is not (solely) due 
to scattered light. Given the long-term activity trend in 
the host sta r (indicating the presence of a magn etic field) 
observed by 
possible exp 



( |2011[|2012[ ), another 



Gomes da Silva et al. 

anation is that the planet induces a spot on 
the stellar surface, which remains synchronized with the 
planet's orbital motion. The fact that the 2007 photom- 
etry does not clearly show a similar strength signal at 
the same period may be related to the magnet ic activity 
cycle detected by Gomes da Silva et al. (2012 1. 

The most convincing case of magnetic star-planet in- 
teractions so far is the Tau Bootis system. A period- 
icity close to the orbital period of the planetary com- 
panio n was detected in M OST photometry of the sys- 
2008[ ), indicating the possible exis- 



tem (Walker et al. 



tence of stellar spots and/or active regions induced by 
Tau Bootis b. Spectropolarimetric observations suggest 
the magnetic polarity reversal may b e at least partly re 



lated to the presen ce of the planet (Donati et al. 



2008 

Fares et al1f2fj08|). GJ 581b (Mp'sini = 0MM Jup ; 
Forveille et al.|[20"ll| ) is less massive than Tau Bootis b 
{M P sin » = '6.9M Jup ; |Butler et al.||1997[ ). However, the 
same can be said for the host stars {M\ = 0.3M Q and 
1.3M for GJ 581 and Tau Bootis, respectively). Fur- 
ther, the b planet is also the largest in the GJ 581 system. 
Thus, star-planet interactions are a viable possibility but 
further monitoring of the GJ 581 system is necessary to 
better understand the nature of the photometric variabil- 
ity observed in the 2009 MOST photometry. 

7. CONCLUSIONS 

We observed GJ 581 with the MOST space telescope 
for a duration of 28 days in 2007 and 11.5 days in 2009. 
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Fig. 10.— The 2007 and 2009 photometry phased individually 
(top and middle panels) and combined (bottom panel) at the or- 
bital period of GJ 581b. The black solid points represent the pho- 
tometry binned every 0.005 orbital phase. The red open squares 
correspond to the data averaged into 0.05 phase bins. 



The 2007 photometry served to determine the level of 
variability of the host star, while the higher-cadence 2009 
photometry was collected with the purpose of searching 
for transits of GJ 581e. We computed 3cr transit windows 
using bo th the Keplerian and circular ephemerides re- 
ported in Forveille et al. (2011 ). To quantify the sensitiv- 
ity of our data to super-Earth transits of various depths, 
we performed Monte Carlo simulations which involved 
injecting simulated transits into the light curve and at- 
tempting to recover them using a box-shaped model. The 
simulated transits have planetary radii between 1 and 4 
J?0 and impact parameter values ranging from 0.0 to 0.9. 

Our transit search returns no signific ant candidate 
This is in agreement with the findings of |Forveille et al 



(2011 1. However, our data allow a full search of the 3<r 
circular orbit transit window, and the 3a Keplerian orbit 
window as well, ensuring that the potential transit of G J 
581e (if sufficiently deep) was not missed by our obser- 
vations. We place lower limits on the planetary density 
of GJ 581e, if it transits. If it docs not, then its mass 
likely lies between 1.84 M ffi and 3.03 M$ (assuming or- 
bital co-planarity and stability over more than a few My r 
for planets b, c, d and e) as found by Mayor et al. ( 2009), 
but we cannot constrain its density. 

We conclude that if GJ581e has an envelope dominated 
by water, helium and/or hydrogen (as appears to be the 
case for GJ 1214b and Kepler- 1 If), we rule out transits 
for most geometric configurations (b < 0.6). We can- 
not rule out transits for a higher density planet. Given 
the significance of this planetary system, we recommend 
that RV measurements continue to be collected in or- 
der to constrain the eccentricity of the orbits more pre- 
cisely. Subsequently, further photometric observations 
should be undertaken to conclusively determine whether 
GJ 581e, if it has a density higher than 2.5-3 g cm~ 3 , 
transits. A positive result would provide much needed 
insight into the bulk and atmospheric composition of this 
low-mass super-Earth. 

Since the MOST photometry spans two orbital cycles 
of GJ 581b, we are able to rule out transits for all plau- 
sible compositions of this planet, but only for low im- 
pact parameter (6 Js 0.4). This leaves a significant range 
of unexplored transiting configurations. We recommend 
the same course of action as for planet e: once the orbital 
parameters of the system are more tightly constrained, 
combining existing GJ 581 photometry with a few hours 
of new strategically acquired photometry should allow a 
definitive conclusion to be drawn regarding the presence 
or absence of GJ 581b transits. 

Transit searches for GJ 581c and GJ 581d cannot be 
adequately carried out with the photometry presented 
in this paper, though we are able to place a very loose 
constraint on the size of GJ 581c, if it transits, in the 
form of an upper limit of Rp — 2.3i?0. 

By combining the 2007 and 2009 photometry, we find 
that GJ 581 is stable to within 1%, indicating that it 
is likely a quiet M dwa rf and supporting evidence of its 
advanced age (> 7 Gyr; Selsis et al. 2007 1 . This increases 
the probability that life may occur and evolve on a planet 
in the habitable zone of the system. We attempt to verify 
the stellar rotatio n period of 94.2 days determined by 
Vogt et al. (2010), but our data set does not cover a 
sufficient portion of a rotation cycle of this length to 
allow any definite conclusions to be drawn. Finally, in 
the 2009 data set we detect a signal with a period of 5.586 
± 0.051 days, which is close to the orbital period of GJ 
581b. A signal at the same period and of similar strength 
is not present in the 2007 photometry. We phase the two 
data sets at the period of the planet, individually and 
combined, and plot them in Figure 10. Though a signal is 
not apparent in the 2007 light curve, one is clearly visible 
in the 2009 and the combined light curves. Additional 
observations are needed to shed more light on the nature 
of this possible phase variation. 
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